Coulomb drag as a measure of trigonal warping in doped graphene 
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I suggest to use the effect of Coulomb drag between two closely positioned graphite monolayers 
(graphene sheets) for experimental measurement of the strength of weak non-linearities of the spec- 
trum in graphene. I consider trigonal warping as a representative mechanism responsible for the 
drag effect. Since graphene is relatively defect-free, I evaluate the drag conductivity in the ballistic 
regime and find that it is proportional to the fourth power of the warping strength. 



The first experimental measurement of conducting 
properties of graphene [ij (an atomically thin crys- 
talline monolayer of graphite) was followed by develop- 
ing graphene-based transistors @, Q , where high concen- 
trations of charge carriers can be induced by applying 
gate voltages. These discoveries have brought a lot of 
attention to the field, which is well-studied theoretically. 
Indeed, a two-dimensional, hexagonal lattice of carbon 
atoms is a usual starting point for most calculations on 
bulk graphite, carbon nanotubes [4|, or fullerenes [5|. In 
fact it's been almost 60 years since the band structure of 
graphene has been first studied 6]. 

The continuous theoretical interest in graphene is due 
to the Dirac-type dispersion relation 0, |I[ leading to a 
number of peculiar properties - from the Berry phase 
of electronic wavefunctions Q to an anomalous Quan- 
tum Hall Effect The chiral nature of charge carri- 
ers in graphene is a consequence of its crystal structure. 
The honeycomb lattice contains two equivalent sublat- 
tices. Nearest-neighbor hopping between A and B sites 
results in formation of two energy bands which intersect 
near the corners of the he xag onal Brillouin zone [3] ■ Two 
inequivalent corners K± [10( define a valley index, speci- 
fying excitations around the Fermi energy (playing a role 
of a pseudospin). Close to the crossing points the spec- 
trum is conical with the "light-speed" v w 10 6 m/s 




FIG. 1: A Coulomb drag sample: two graphene sheets are 
separated by an insulating layer of the thickness d. The two 
gates at the top and bottom of the device can be used to 
independently control carrier concentrations in the two sheets. 

The true microscopic Hamiltonian in graphene con- 
tains several small corrections to the Dirac spectrum. For 
example, using a second order k • p equation one can 
derive a quadratic term in the effective low-energy Hamil- 
tonian of graphene that violates the isotropy of the Dirac 
spectrum and causes trigonal warping [9(. When gate 



voltage is applied, resulting in non-zero Fermi energy, the 
electronic Fermi line deviates from a perfect circle: the 
p — > — p symmetry of the Fermi surface is broken within 
each valley. Although weak, such distortions of the Fermi 
surface lead to observable effects: trigonal warping sup- 
presses antilocalization [lit ], which one would otherwise 
expect due to the absence of backscattering [3, [l3| . 

Trigonal warping also breaks another, more subtle 
symmetry in the problem - the symmetry between exci- 
tations above and below the finite Fermi energy Ep = 
hvy<irn e (n e being the electron density). Such asym- 
metry can be detected in transport measurements. A 
particular "tell-tale" experiment that crucially depends 



on its presence is Coulomb drag [1J, [15|, 

lid na us 

[Til , lifjj ]. Coulomb drag measurements are performed 
on two closely positioned (but electrically isolated) lay- 
ers. A current I a driven through one of the layers (the 
"active" layer) induces a voltage drop V p = pr>Ia in 
the other ("passive") layer. The voltage appears due 
to inter-layer electron-electron interaction that creates a 
frictional force "dragging" electrons in the passive layer. 
As a purely interaction effect, Coulomb drag has be- 
come a sensitive tool for experimental studies of electron- 
electron interaction in many problems of contemporary 
condensed matter physics. It has been used in search for 
Bose condensation of interlayer excitons 2l|, a metal- 
insulator transition in two-dimensional layers [22| . and 
Wigner cristallization in quantum wires [231 ] . 

Other effects contributing to non-linearity of the spec- 
trum in graphene will also result in non-zero contribution 
to the Coulomb drag. These include: the quadratic cor- 
rection to quasiparticle spectrum due to next-neighbor 
hopping Coulomb scatterers resulting in energy- 
dependent scattering time r c oc \Ep\ 12411 : logarithmic 
corrections to quasiparticle spectru m [25l| : and interfer- 
ence corrections to scattering time [26J. In this Letter I 
consider trigonal warping as a representative mechanism 
of the Coulomb drag in graphene, leaving the discussion 
of the role of other mechanisms for a subsequent publi- 
cation [27| ■ I argue that this effect can be distinguished 
from other drag mechanisms by its dependence on inter- 
layer separation d and Fermi momenta (or gate voltages). 

I envision a following set-up (see Fig.QJ. Two graphene 
sheets are positioned parallel to each other and are sep- 
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arated by an insulating material about 50 nm thick (e.g. 
using the technique recently developed in Ref. [28h . Two 
gates (at the top and bottom of the device) can be used to 
independently control carrier concentrations in the two 
layers. I assume that sufficient gate voltage is applied, so 
that Fermi energies in each layer are positive Ep P ' > 
and represent the largest energy scale in the problem (in- 
dices a and p denote the active and passive layers). Since 
graphene is relatively defect-free, the elastic scattering 
rate r _1 is assumed to be much smaller than temper- 
ature. At the same time, the mean-free path i is the 
longest length scale. The assumed hierarchy of energy 
and length scales (here Xf is the Fermi wavelength) 

hr- 1 « T « 4 Q ' P) - Z (a ' p) » d » \f p) , (1) 

ensures that the the device is in the ballistic regime (29[ 
(for Coulomb drag in the ballistic regime see Ref. flil ) . 

At lower temperatures (T < hr^ 1 ) electron motion 
becomes diffusive. In that case, scattering off atomically 
sharp disorder becomes important for it breaks the pseu- 
dospin symmetry and drastically affects two-particle cor- 
relation functions [13] • An analysis of the Coulomb drag 
in the diffusive regime will be considered elsewhere 27 1. 

In the weak coupling regime [HI, [l^, [2(| the drag co- 
efficient is proportional to the drag conductivity « 
ujy/(a a a p ) (<J( a .p) being the Drude conductances of the 
two layers). The latter is typically calculated using the 
expresssion 
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where S is the sample area, T> ap is the screened inter- 
layer interaction, and T is the non-linear susceptibility 
(or rectification function) that relates a scalar potential 
V(u>, q) to the dc current it creates in quadratic response 
j dc = r(w,q)|l/(w, q)| 2 . Below I re-derive Eq. © for 



graphene in the ballistic regime under the above assump- 
tions and show, that for Dirac particles the drag vanishes. 

When trigonal warping is taken into account, I find 
that the drag conductivity is proportional to the fourth 
power of the parameter W that describes the strength of 
the warping correction to the Dirac spectrum: 
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Here Kt a ,p) = e 2 kp' p ^ /v are the Thomas-Fermi momenta. 
Eq. ([3]) is the main result of this communication. 

The low-energy single-particle Hamiltonian 0, 0, [13 
in graphene can be written in the basis of 4-component 
Bloch functions $ = (<j)A.K + , 4>b.k + , (t>B k_ , <t>A.K_ ) (I 
employ notations introduced in Ref. 12; 4>a,k + is the 
electronic amplitude on sublattice A and valley K+) as 



Ho = wEp - W 



(f 2 



{Px - Pi) ~ ZVyPxPy 

with the weak quadratic term causing trigonal warping. 
Here Pauli matrices &i act in the sublattice space (A, B). 
The "isospin" E is defined as direct products of Pauli 
matrices a (acting in the sublattice space) and II (acting 
in the valley space K±): ^ x i y ) = H Z ® &x(y)- 

In the basis of plane waves H.q is a 4 x 4 matrix that can 
be diagonalized by a unitary transformation i? _1 7io-R = 
diag[£' Q ^(p)]. The resulting eigenvalues are 



(4) 



E a ,£ — avps^: = y 1 — 2£pWv~ 1 cos3ip p +p 2 W 2 v 2 , 

(5) 

where £ = ±1 denotes the two valleys, a = ±1 is the chi- 
rality and distinguishes between the conductance (a = 1) 
and valence (a = — l)bands, and <p p is an angle between 
the momentum p and the x-axis (tan</? p = p y /p x )- 

The electron field operator can be written in the basis 
of the eigenstates as (hereafter I use the units with H = 1) 
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where e a ^ — as^e 2liPp /[e 3llf,p — (£pW/v)]. Then the form 
of the electron density operator 

p,p',a,a' ,^ 

differs from the usual one by the presence of vertices 

A££ = (l + e«^)/2. (8) 



The vertices A k ' k , indicate the asymmetry of quasi- 
particle scattering in graphene. In particular, the su- 
pression of backscattering [9] follows from the fact that 
in the absence of trigonal warping Xp'" p (W — 0) = 0. 
Similarly, Xp'p a (W = 0) = 0, indicating that inter-band 
transitions at the same wave-vector are also supressed. 

The single-particle Green's function in the original ba- 
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sis of Bloch functions is a 4 x 4 matrix [12j It can be 
brought to a diagonal form by the rotation R. In any 
closed loop such operation would cause the appearance 
of the vertex factors (jHJ). For example, the polarization 
operator (the density-density response function) is 



P (w,q) = -2» 



d 2 k 

(27T) 2 



de 
2^ 



(9) 



k,k+q 



G a< £ (e, k) G Q '^(e + w, k + q) 



where the overall factor of 2 follows from the spin de- 
generacy. Here G Q (e, k) is the Green's function for the a 
band. For Ep > and T <C Ep one can use the standard 
non-relativistic Green's functions (30j . 

Since the vertex functions are frequency- independent, 
the usual reasoning [3l[ leads to the following expres- 
sion for the retarded polarization operator in the ballistic 
regime (f[E a (k)] is the Fermi function) 
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where r/ — > +0. Note, that under assumptions that the 
Fermi level is in the conduction band and T <C Ep, the 
term with a = a' = — 1 (the valence band contribution) 
vanishes due to the Fermi functions. 

For the purposes of describing screened inter-layer in- 
teraction in the Coulomb drag problem, I am only inter- 
ested in momenta q smaller than inverse inter-layer dis- 
tance. Under the assumption (JXJ) , q < (1/d) <C kp. Thus, 
the inter-band contribution to Eq. (|10[) is suppressed by 
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a factor of q 2 /kp and the polarization operator is dom- 
inated by the conduction band. Due to the Fermi func- 
tions in Eq. (|10[) the momentum integral is dominated by 
the region k ~ kp. Then the leading contribution to the 
static polarization operator needed in the Coulomb drag 
problem in the ballistic regime [3] is 



0) = 2k F /irv = 4i/. 



(11) 



Here v is the density of states at the Fermi level (spin and 
valley degeneracy is taken into account). Consequently, 
the screened interlayer interaction is the same as in the 
case of the uniform two-dimensional electron gas [l9| : 



Dab = Tre 2 q/(K a K p smhqd). 



(12) 



To derive the expression for the drag conductivity @ 
one starts with the general expression for electric cur- 
rent in the passive layer in terms of the Keldysh Green's 
function. In graphene the current vertex (in the original 
basis of Bloch functions) is J = 2ewS. Diagonalizing the 
Green's function by a unitary transformation the current 
in the passive layer takes the form 
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In a system of free electrons in equilibrium the current 
(fT3|) is equal to zero. Perturbing the Green's function 
by a potential V(co) one finds the following expression 
for the current [more precisely, the contribution relevant 
for the drag problem; here I use a short-hand notation 
for the spatial coordinates in the argument of Green's 
functions - G£ £ (e; 31) = G^(e; r 3 , n)] 
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In the ballistic regime [19( Green's functions in Eq. (TTj 
are averaged over the disorder independently of each 
other. Then one can perform the Fourier transform and 
use the approximation (r is the elastic scattering rate): 

G^(e; k)G^(e; k) » 2nrS(e - £7 a , £ (k)) , (15a) 



ImG^(qk) w-7r*(e-£? a , c (k)). (15b) 
In this case 

f[e] - f[e- iv] -» f[E a ,t(k)] - f[E a ,t(k - q)]. 



It is then clear that if the Fermi level is in the conduc- 
tion band then the two Fermi functions for the valence 
bands cancel each other. Inter-band processes are sup- 
pressed similarly to the case of the polarization operator 
and will be neglected hereafter. Thus, in the ballistic 
regime, only particles in the conduction band contribute 
to the current in the passive layer, as one would intu- 
itively assume. The question of whether this statement 
remains true when off-diagonal disorder is taken into ac- 
count, i.e. in the diffusive regime, will be discussed in a 
subsequent publication [27| . The situation in the active 
layer is similar, in fact both layers are described by the 
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same non-linear susceptibility T. Therefore, the general 
expression for the drag conductivity ([2]) remains valid 
(again, with only particles in the conduction band con- 
tributiong). What remains to be done inorder to obtain 
the result ^ is to evaluate the non-linear susceptibility. 

Under my assumptions the non-linear susceptibility for 
the conduction band in the ballistic regime is 



UJ 



r(w,q) = - PfqM - 7-q(-w)] 



where the triangular vertex 7 is given by 
d 2 k 



1.1 

k,k+q 



(16) 



(17) 



x [lmG^(e + cj; k + q)] (e; k)] x G^ (e; k) . 

Here j 1 is the diagonal matrix element of the curent ver- 
tex J rotated (|13[) to the basis of the eigenfunctions: 
ji = 2evrik. Note how the current vertex for the con- 
duction band recovers its usual momentum dependence! 

Consider now the Coulomb drag for Dirac particles. 
Setting W = in Eq. ©, one finds [9 = Z(k,q)] 
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The result is even under the simultaneous change of sign 
of u) and q: 7 q (w) = 7- q (— w). Therefore r(w, q) = 0. 
There is no drag effect in a system with linear spectrum. 

When the deviation from linearity in Eq. (Tj[]) is taken 
into account, the contribution of the two valleys to 
Eq. I|18p is no longer identical and a non-zero result ap- 
pears only in the second order in W: 

r(cj,q) = — 4ei , qujw 2 £v~ 3 sin</? q cos 3<,Sq#(v<7 — uj), (19) 

where the mean-free path is defined as I = 2vt. Using 
Eqs. (TO]) and JT5J in Eq. (©, I find the final result ©. 

To summarize, I have considered the Coulomb drag be- 
tween two closely positioned graphene sheets. For strictly 
linear Dirac-type dispersion, I find that the drag vanishes, 
in agreement with the traditional interpretation of the ef- 
fect as a manifestation of asymmetry between elementary 
excitations above and below the Fermi level. As a rep- 
resentative mechanism of such asymmetry in graphene 
I consider trigonal warping and find the drag coefficient 
proportional to the fourth power of the strength W of 
the warping term [32j. The obtained result should be 
distinguishable from the drag due to other non-linear 
contributions 0, 24, [SB], [26| to graphene spectrum by 



its dependence on inter-layer separation and Fermi mo- 
mentum. In my opinion the Coulomb drag is an ideal 
tool for experimental studies of spectrum non-linearities 
in graphene. 
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